Abstract 

A numerical scaling analysis is used to show that Nagaoka's ferromagnetic 

state in two-dimensional Hubbard model with one hole is supersede by an 

antiferromagnetic (AF) state with a discontinuous jump in the total spin due 

to the AF coupling as the Hubbard U is made finite. The same applies to 

t::J" ' the two-hole system, which has a spiral spin structure. We can show, via 

^^ ■ the scaling, that the crossover to an AF state is a precursor of a pathological 

^ ■ coalescence of states having the minimum spin and Nagaoka's state at f7 = oo 

f-^ ■ in the thermodynamic limit. 

PACS numbers: 75.10.Lp, 71.27.-ha, 75.30.Ds 
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To identify a stable itinerant ferromagnetism is one of the longest- standing yet unresolved 
problems in the physics of strongly correlated electron systems. Nagaoka's theorem is the 
first rigorous result, which served as a foundation of subsequent studies. [|1],0,0] Nagaoka's 
problem still holds a fascination as a remarkable manifestation of the correlation effect in the 
doped Mott's insulator, where a key concept is interference between the motion of a hole and 
its surrounding spin configuration. 0] A recurring question on Nagaoka's ferromagnetism, 
however, is the singular condition (the half-filled band doped with a single hole with the 
Hubbard U = oo): Is it possible to extend the ferromagnetism away from the extreme 
condition down to finite number of holes and/or finite coupling strengths? 

There have in fact been many attempts at looking into the possibility of a finite fer- 
romagnetic region in the vicinity of Nagaoka's limit in the phase diagram. Nagaoka has 
already shown that the spin-wave excitation in Nagaoka's ferromagnet is singular in that 
the spin stiffness vanishes like the inverse of the system size, A^. If one concentrates on the 
stabihty against one-spin flip (ISF) modes having 5* = ^max ~ 1, however, the ferromag- 
netic phase seems to extend over a flnite area on the parameter space of the density of holes 
{6 = Ni^/N) and the antiferromagnetic (AF) coupling (J = At'^/U), [§|,§|,0,§|-0 but these vari- 
ational (usually Gutzwiller-projected RPA) results should be considered as an upper limit 
for the ferromagnetic boundary. On the other hand, a recent high-temperature expansion 
result by Putikka et al. shows that the true ground-state is signiflcantly lower in energy 
than the fully-polarized state for any hole density. 0,0 They also argue the possibility 



of a flnite 'ferrimagnetic' region, which is identifled from the uniform spin susceptibility 
diverging like the Curie law. 

An exact way to explore the problem is the numerical diagonalization of flnite systems. A 
serious point found in early studies is that Nagaoka's state is sensitive to both the boundary 



condition and the number of holes, [n,|T3 where singlet states can supersede Nagaoka's state 



even aX, U = oo. These results have motivated subsequent investigations. OJT^JTOJT^ As 



for the instability of Nagaoka's state against J, loffe and Larkin |]TB[ have conjectured that 



a crossover to an AF phase with a jump in S can possibly occur for increasing J . Nature of 



the lower-spin states that compete with Nagaoka's, however, is not fully understood yet. 

The purpose of the present paper is (i)to give the low-lying spectrum including both 
charge and spin excitations, which will serve as a foundation (ii)to obtain the scaling prop- 
erties of the energy gap and level crossings to actually identify how Nagaoka's state is 
superseded for the 2D square lattice. In the formalism for (i), the heavy mass spin- wave, 
which causes the coalescence of states comes from an interference between the charge (Bloch) 
momentum and the spin momentum. In (ii) we move on to calculate the low-lying states 
for every possible total spin (0 or 1/2 < S < Smax) to show that the first state to take over 
Nagaoka's as U is decreased from infinity is indeed an AF state rather than a partially polar- 
ized one, which indicates a discontinuous jump in the total spin. The same sudden crossover 
occurs for the S = ground state at U = oo that appears for one hole with the anti-periodic 
boundary condition (APBC) or for two holes with the periodic boundary condition (PBC), 
which we identify as a spiral spin structure with a wave-length comparable to the system 
size. The crossover by an AF state is related, via the scaling, to a pathological degeneracy 
at f/ = oo, where even the lowest-spin states become degenerate with the Nagaoka's (or 
spiral) state in the thermodynamic limit (A^ -^ oo). 

In order to grasp the subtleness of Nagaoka's state, we start with introducing the follow- 
ing formulation. [|l9l Consider the laxge-U Hubbard model on a (i-dimensional hyper-cubic 
lattice with PBC. An eigenfunction for N^. = N — 1 electrons with one hole can be written 



as 



10) = ^exp(ifc-a.h)E/(^i'---'^^c)^k<..---^L. <...)' (0-1) 

where x^ is the position of a hole and Cxa- creates an electron with spin aj at Xj. The 
coefficient, /(o"i, ■ ■ ■ , cr^J, can be regarded as a spin wavefunction around the hole. Note that 
we do not sum over the position of each electron: A trick in writing the above formula is that 
the sequence (a^i, . . . , x^J (where no two positions coincide to exclude double occupancies) 
of the multiplication of creation operators is fixed in such a way that the creation at the 
same coordinate relative to aih enters in the same position for each value of a^h- In this 



convention / no longer contains a^h reflecting the translational symmetry, while the factor 
exp (ifc ■ a^h) takes care of Bloch's theorem. 

Consider the strong-coupling model, which is derived from the Hubbard model in the 
second order in l/?7-expansion. Operating to the above basis we end up with a spin Hamil- 
tonian, 

d 

^spin = -t "^{exp {ik ■Ua)Ta + exp {-ik ■ Ua) T^^] 
+ J E (5, ■ 5, - \) 

-jJ2 ^^P (^^ -Ua + ik- Ua') TaTa> 



a, a' 



X 



[riua+u^, ^riu^, I + nu^+u^, inu^, t " ^Ua+u^, ^u^, - Su^+u^, Su^^ ) , (0.2) 



which corresponds to the spin part of the t — J model (including three-site terms). Here 
the operation, 7^(a = a;, y, ...), is roughly a translation of the spin configuration along a-axis 
except for a spin on XY^ + Ua-, which is shifted onto a^h ^ '"a with Ua being a nearest-neighbor 
vector. The last line represents the three-site term, where the summation involves the 
neighbor {ua) of the hole site and the third site {ua') which is the neighbor of Ua- 

The expression in the first term of Eq. ( P^ ) suggests the following picture. Since 7^ is 
almost a translation, it gives a phase which is roughly the momentum of the spin configu- 
ration. To be more precise, we can introduce, in place of /(o"i, ■ ■ ■ , cr^r^) where spins are not 
allowed to occupy the hole site {i.e. Xj ^ 0), another basis, g{yi, ■ ■ ■ , Vm)^ ^^ M-spin-fiip 
states on the square lattice. Here t/^, which may be 0, denotes the position of j-th down 
spin. We can then readily construct an eigenstate of the total momentum, Q. As a penalty 
for using the extended basis, the operation Ta must be redefined: it gives zero if any y, 
equals to in g. In this representation, the problem becomes a scattering of a spin at the 
hole site. Because the scattering is a higher-order process of 0{1/N) compared with the 
translation, we can obtain the asymptotic expansion in 1/A^ for the ISF heavy-mass mode 
in 2D. The solution exhibits a heavy mass (the band width ~ 4'7rt/A^ln A^), which coincides 
with the result by Barbieri et al. [^ 



The advantage here is that the present formahsm enables us to introduce a quasi- 
momentum Q defined by Ta\(j)) = exp{iQ ■ Ua)\(f)). Then the asymptotic resuh indicates 
that Q is nearly equal to the momentum, Q, of the incoming state (with \Q — Q\'^ ~ 
0(l/A^lnA^) for spin-wave excitations). For a solution with Q{k), the dispersion is given 
by e{k) = —2tY,'!i=i^'^^iQaik) — ka) for each k. In a finite system we have discrete Ap- 
points (= 0, fci, . . .) having single-particle excitation energy e{= 0, ei, ...). From this we can 
predict that the ISF excitation from Nagaoka's state should comprise a series of bands, 
each of which consists of an incoming wave, Q = k + ki and has a quasi-momentum 
Q = k + ki + Ak{i = 0,...), which are separated by the single-particle energy, £,, with 
the band widths vanishing like 0(t/A^lnA^). Note that Ei itself vanishes like 0{t/N), which 
roughly gives the gap between heavy-mass bands. 

This picture is in fact confirmed from the numerical diagonalization of a 290-site system 
in Fig. |l[ Interestingly the above picture also holds for two-spin-fiip (2SF) excitations, where 
the width of the lowest 2SF continuum as well as the width of the lowest ISF dispersion fit 
to a size dependence of Ai5i7(oo), A25_f(oo) ~ t/N\nN ^ £1 ~ 0(t/N) within a finite-size 
correction (inset of Fig. |I|). The energy of 2SF mode is shown to be approximately twice 
that of ISF, so that the spin-waves are nearly free. This is natural because the interaction 
between spin-waves is a higher-order process than the scattering of a single spin at the hole 
site. 

We now turn to finite U for a fixed boundary condition in finite systems to probe level 
crossings between Nagaoka's state and other states having various total spin. So far Na- 
gaoka's state is suggested to change to partially polarized ones, then finally to a lowest-spin 
state after several level crossings in 2D from numerical studies on 10- and 16-site systems. 
EI|,|22|,g|,|2§ 



However, if we increase the size to a 20-site system here, the strong-coupling model {t-J 
model with or without the three-site terms) shows the absence of the intermediate partial 
polarization in the ground state, (Fig. ^ where the range of J for the partially polarized 
region decreases systematically to zero (Fig. §). This suggests that an abrupt crossover from 



the fully spin-polarized state to an unpolarized state occurs for larger systems. 

We can relate this feature with the spectrum in the strong-coupling limit via the scaling 
property of the level crossings. We have previously reported from the value of f/c for the 



level crossings that a 2SF mode takes over Nagaoka's state before a ISF mode does so. ||T^ 
If we plot here the crossing point, U^^ , between Nagaoka's state and the AF state as well, 
we see that this is the crossing that comes downwards most rapidly for increasing J (inset 
of Fig. 1^). Thus we obtain a lower bound for U^^ as U^^ > —NlnN. The asymptotic form 
for the energy of the AF state is expected to be 

Eaf = -4t + Aaf - aN— , (0.3) 

where Aaf is the energy gap between the AF state above Nagaoka's state at f/ = oo and 
the last term with a ~ 0(1) is the AF exchange interaction. Hence we end up with an upper 
bound for Aaf, via Eaf{U^^) = -ENagaoka = -4t, as 

Aaf < ant/ \nN , (0.4) 

which vanishes in the thermodynamic limit. Therefor the huge degeneracy of the ground 
state in Nagaoka's limit is seen to include even the AF state. 

Although it is premature to discuss whether a spin-polaron picture or string picture HJ25[ 



applies here, we have confirmed that the lowest 5* = 1/2 state has the energy that varies 
approximately linearly with J for J^^ < J < O.lt and S{Q) that remains to be peaked at 
(7r,7r), which implies that the AF state experiences no level crossing prior to the crossing 



with Nagaoka's for systems up to 26-sites. (See also p6| .) 

As far as finite systems are concerned, we find that the number of holes is crucial while 
the concentration of holes is ill-defined in the Hubbard model at f/ ~ oo. As already noted 
by Riera and Young, [|13| there are similarities among the systems having the same number 
of holes. Thus we can look for a scaling property in the series of finite samples with the same 
number of holes. For two holes, we can look at the level crossing between the spiral-spin 
ground state for large U characterized by S{Q) having four peaks around the F-point and 

6 



the AF state. If we fit the crossover point from Nagaoka's state, J^^ (= 4£^ /U^^) against 5 
with a power-low, we have Jl^°^° oc 5^'^^, J^'^°^'^ oc ^^■^^. Thus for J ^ Nagaoka's (or spiral) 
state is only realized for an infinitesimal doping, if this scaling persists to 5 ^ 0. We notice 
in Fig. ^ that two scaling curves, although having similar exponents, differ considerably, 
which endorses the difficulty in defining the hole concentration. Since both 7^^*^°'^ and J'^^'^^^ 
are tangent to the 5 axis, a convex functional form for other number of holes would be 
required to have a finite region for magnetism. If we compare our scaling relation with the 
result (dashed line in Fig. |^) from the high-temperature expansion by Putikka et al, the 
width for the partial polarization for one- and two-holes vanishes in what they call is the 
ferrimagnetic region. As far as finite systems with A^ ~ 20 are concerned, our result thus 
indicates the absence of magnetism in this region. 

However, since a finite number of holes only correspond to an infinitesimal hole doping 
for A^ — * oo, the true problem of the ground state for finite concentrations of holes in the 
thermodynamic limit remains an open question. 

We wish to thank Prof. M. Ogata for illuminating discussions. We are also grateful for 
Prof. Y. Nagaoka in the early stage of the present work. 
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FIGURES 
FIG. 1. Spin- wave excitations from Nagaoka's state (solid circle) in a 290-site 2D square lat- 
tice with U = oo. Open circles represent one-spin-flip (ISF) excitations, while hatched regions 
represent continuum of two-spin-flip excitations. The arrow denotes the first single-particle ex- 
citation energy, Agmgie- The inset shows the size dependence of the band width, Aisf, of the 
lowest ISF mode (solid line), the width, A2SF1 of the lowest 2SF continuum (broken line), and the 
bottom-to-bottom spacing of the first two 2SF continua, IF2SF (chain line). We also indicate the 
dependence, 47rt/A^lnA^ (lower dotted line), and the single-particle energy Asingic('^ 1/^) (upper 
dotted line). 

FIG. 2. Level crossings in a 20-site t-J model in the small J regime. The ground state changes 
from Nagaoka's state to a singlet state at J^^ = 0.0413t. The inset shows the scaling of various 
C/c's for t-J models with 3-site terms defined in the text. The dotted line represents the analytic 
evaluation, U^^^ = {4:t/7r)NlnN, while the dashed, chain and solid lines are a guide to the eye. 

FIG. 3. Scaling for the crossover point, Jc, for the one-hole doped case between Nagaoka's state 
and the antiferromagnetic (AF) state, and for the two-hole doped case between the spiral-spin state 
and the AF one. Bars represent not the error bar but the range over which partially polarized 
states are realized. A broken line is the boundary for the 'ferrimagnetic' region obtained by Putikka 
et al. (Ref. 0,|ll|). 
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